Incipient nodal pairing in planar d-wave superconductors 
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The possibility of a second pairing transition d —> d + is {d + id') in planar d-wave supercon- 
ductors which occurs in the absence of external magnetic field, magnetic impurities or boundaries 
is established in the framework of the non-perturbative phenomenon of dynamical chiral symmetry 
breaking in the system of 2-f 1-dimensional Dirac-like nodal quasiparticles. We determine the critical 
exponents and quasiparticle spectral functions that characterize the corresponding quantum critical 
behavior and discuss some of its potentially observable spectral and transport features. 
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The recent ARPES ||I|] and optical conductivity 
studies of the superconducting high-T^ compound 
Bi2Sr2CaCu20s+s have chaUenged the physical picture 
of what was once thought to be an essentially mun- 
dane BCS-like pairing state with only weakly interacting 
quasiparticles near the nodes of the d^2_y2 order param- 
eter. Most notably, in a striking contrast with the well- 
known theoretical prediction of the cubic in temperature 
quasiparticle damping rate obtained in the conventional 
d-wave picture [p| , both of the above probes plO] revealed 
a much weaker, approximately linear, temperature de- 
pendence. 

Whether the above findings indicated a behavior 
generic for all cuprates or not, their observation has 
prompted the idea of a quantum-critical behavior associ- 
ated with a zero temperature quantum phase transition. 
Should such a transition occur in a sufficient proximity 
to the physical slice (the standard " temperature-doping" 
plane) of the, conceivably, multi-dimensional phase di- 
agram of the cuprates, the fluctuations of a new order 
parameter could provide for a strong enhancement of 
quasiparticle scattering resulting in the sought-after lin- 
ear quasiparticle damping rate. 

Recently, M. Vojta, Y. Zang, and S. Sachdev pro- 
posed a phenomenological theory of the nodal quasipar- 
ticles coupled to the fluctuations of a secondary (is or 
idxy) pairing order parameter. By applying a one-loop e- 
expansion about D = 3 + I space-time dimensions, they 
found an evidence of a quantum critical point residing 
inside the superconducting phase Q. 

In the present Letter, we develop a direct 2-1-1- 
dimensional approach that will allow us to establish the 
existence of the critical point in question beyond pertur- 
bation theory and make quantitative predictions for the 
corresponding quantum critical behavior. 

In order to facilitate the following discussion we make 
use of the rotationally-covariant Nambu representation 
where the d-wave quasiparticles are described in terms of 
the two species (i = 1,2), each having N — 2 spin compo- 
nents labeled by a, of the two-component Nambu spinors 
V'ia = (cQ;(ki), eQ,/3cL(— ki)) which are constructed from 
the quasiparticle states near the pairs of the opposite 
nodes located at ±ki = ±{kp,±kp)/\/2. 

Upon rotating the axes of the coordinate system by 



7r/4, we cast the linear in momentum kinetic energy of 
the nodal quasiparticles in the form of the (anisotropic) 
Dirac Lagrangian 
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where vp and wa are related to the momentum deriva- 
tives of the normal state quasiparticle dispersion and the 
parent (i^2_j^2 -symmetrical gap function, respectively. In 
Eq.(l), we use the irreducible 2x2 representation for the 
7-matrices in 2 -I- 1 dimensions 7^ = {f72,iiyi,i(J^) and 
the notation -i/^j = Wilo- Also, in order to make Eq.(l) 
more symmetrical we perform a rotation in the Nambu 
space on one of the two fcrmion species: 'ip2a — >■ '02a = 

((Ti -f Cr3.)lp2a/\/^- 

Next, we combine the spinors i/jic and "020: into one 
four-component Dirac fermion '^a. = i'4'ia,'4'2a) and in- 
troduce a (reducible) 4-dimensional representation for 
the 7-matrices r^ = 7^ (g) CT3 where the second factor 
in the tensor product acts in the subspace of the two 
species. 

In the phenomenological secondary pairing scenario, 
the nodal Dirac fcrmions couple to a real boson field (f> 
corresponding to the imaginary part of the total gap func- 
tion. Its own critical dynamics is governed by a generic 
0*-theory 
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At T = the theory (0) undergoes the conventional D = 
3 Ising transition driven by the quartic self-interaction 
of the field that results in spontaneous breaking of the 
reflection {Z2) symmetry (f) -^ —(j) as the parameter m?' 
is tuned into its critical value m^. 

Alternatively, the Z2 symmetry breaking can also be 
driven by a linear coupling between (f> and the nodal 
fermions 
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provided that both (mj and (0) remain invariant under 
the discrete chiral transformation 



^c^Ts^fa, 
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where Fg = — 1 (g) (72 anticommutes with ah F^. 

By choosing A = 1 1 in (Bh, one couples (f> to a. 
Z2-(and, accordingly, time reversal-) odd fermion mass 
operator. In terms of the original nodal quasiparticles, 
the latter reads as ipicr2'4'i + W2'^2'4'2 which appears to 
coincide with the order parameter of the is pairing be- 
tween the opposite nodes, while the chiral reflection (4) 
corresponds to a permutation of the two fermion species: 

V'la ^ '02a- 

With all the three velocities in (1) and (2) set equal 
{vp — v/\ = c = 1 ) the theory described by the La- 
grangian C^ + C^ + Cg becomes manifestly Lorentz- and 
Z2-invariant (we will return to the spatially anisotropic 
case later). 

In fact, the above Lagrangian can be readily recog- 
nized as that of the Higgs- Yukawa (HY) model where 
the phenomenon of spontaneous chiral symmetry break- 
ing (CSB) has long been discussed as an alternative to 
the Higgs mechanism [^. Unlike its 1-1-1- and 3-1-1- 
dimensional counterparts which appear to be asymptot- 
ically free and trivial (Gaussian) , respectively, the 2-1-1 
dimensional HY model possesses a non-Gaussian infrared 
fixed point at a finite coupling Qc where CSB occurs 
via spontaneous generation of a non-zero fermion mass 
M = 5(0). 

The above correspondence allows one to identify the 
CSB phenomenon with the conjectured second pairing 
transition (d -^ d ^- is) below which the nodal quasipar- 
ticles become fully gapped. Albeit being absent in any 
finite order of perturbation theory, spontaneous CSB oc- 
curs at the level of the non-perturbative mean field " gap 
equation" for the fermion mass 
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which accounts for the tadpole ("Hartree") contribution 
to the fermion propagator |q|. 

At T = Eq.(^ yields a critical couphng g1 = 
TTTn? /Q,N as a function of the ultraviolet cutoff Vt set by 
the amplitude of the parent d-wave gap. At strong cou- 
pling {g > 5c) 7 the chiral symmetry (4) is spontaneously 
broken, and the fermion mass scales alongside with the 
order parameter Af(g,0) cx< >cx {g — gc)^ ■ 

At finite T the chiral symmetry gets restored above 
a transition line in the g — T plane which terminates 
at the quantum critical point. In the symmetry bro- 
ken (strong coupling, low-T) phase, the mean field equa- 
tion (@) yields the fermion mass M{g,T) = M{g,0) + 
2rin[(l -I- i/l-4exp(-M(g,0)/T))/2] which vanishes 
along a critical line, resulting in the non-BCS-like rela- 
tion M{g,0)/T* = 21n2 between the maximum gap and 
critical temperature T* of the nodal pairing. 



The quantum critical point associated with the CSB 
transition manifests itself via anomalous operator dimen- 
sions that can be found from the solution of the cou- 
pled Dyson equations for the fermion propagator (^^) — 
{Zp^+Ap)~^ with all the tadpoles absorbed into the defi- 
nition of the bare fermion mass (hereafter, -^ = F^p^ and 
p2 ^pt^pp): 



iZp-l)^+Ap = 

d^k Zp_ki^- fi)+Ap_k 
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and the fermion polarization operator 
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which modifies the boson propagator {4>^) = (fc^ -|- 
TO^ -I- Hfe)~^. In particular, the (ultraviolet divergent) 
momentum-independent part of (|^) contributes to the 
renormalized boson mass Sm? = rri^ + Hq which vanishes 
at the critical point g = gc, consistent with Eq.(||). 

In the lowest 1/A^ order, Zp = I, Ap ~ M, Ap^k — 9, 
and the polarization operator (R) assumes the form 



Hfe = Ho + ^^ tan-i 
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In the critical regime {6m = M = 0), the infrared be- 
havior of the boson propagator is governed solely by the 
fermion polarization H^ — Hq ex \J —k'^ (hereafter, ex fc), 
while the bare fc^ term can be completely neglected. This 
implies that at the CSB fixed point the mean field dimen- 
sion of changes from its bare value 1/2 to 1, thereby 
rendering the kinetic and quartic terms in Eq. (0) totally 
irrelevant in the renormalization group sense. 

We mention, in passing, that besides indicating the ex- 
istence of a fixed point at a critical coupling 17^ oc e = 
4 — 13, the one-loop e-expansion used in Ref. jj] also pre- 
dicts that the (formally relevant for any e > 0) quar- 
tic term in (y|) scales towards strong coupling as well: 
A — > Ac oc e, which appears to be an artifact of the cor- 
responding one-loop renormalization group equations. 

With the only relevant term m'^cjp in Eq.(0) remain- 
ing, the critical behavior at the CSB infrared fixed point 
turns out to be identical to the ultraviolet asymptoti- 
cal regime of the 2 -I- 1-dimensional Gross-Neveu model 
which is renormalizable in the 1/iV-expansion, N being 
the number of Dirac fermion species B . 

Namely, the mean field (or N — 00) scaling behav- 
ior of the boson propagator < (j)(f> >cx l/k gives rise to 
the logarithmic divergence of the momentum integrals 
in Eqs.(||,^ describing the first order l/N corrections 



to the fermion and boson wave functions as well as the 
gap and vertex functions. By differentiating these log- 
arithmic corrections with respect to the external mo- 
menta, one arrives at the usual renormalization group 
equations whose solution yields the infrared asyniptotics 
of the renormalization factors 
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from which the dimensions of the fermion [^] and boson 
[(/)] fields can be readily read off. Notably, the anomalous 
dimension of the fermion operator [^] — 1 w l/Svr^iV 
remains rather small even for N — 2. 

As regards the Lorentz-non-invariance of the bare 
fermion Lagrangian (m), in the quantum-critical regime 
the momentum integrals in Eqs.(@,CT are dominated by 
the momenta parallel to the external ones, and, there- 
fore, the anisotropic velocity factors can be scaled away 
without affecting the above power counting. 

The critical exponents characterizing the CSB transi- 
tion satisfy the hyperscaling relations which allow one 
to express all of them in terms of only two indepen- 
dent ones, e.g., the anomalous dimension exponent ry = 
2 - D + 2[0] w 1 - (16/37r2iV) and the exponent ly = 
/3/[0] w 1 + (8/37r^A^) controlling the correlation length: 
1/L^ — 6w? /{g^N) (X \g — gcV ■ The latter remains finite 
in the symmetry broken (ordered) phase, as breaking dis- 
crete chiral symmetry does not result in the appearance 
of a Goldstone mode below T*. 

The fact that the critical exponents demonstrate an ex- 
plicit A'^-dependence implies that the universality class of 
the CSB transition is different for different A^, the D = i 
Ising transition {v = 0.63, r; = 0.025) being recovered 
only in the limit A^ — > 0. 

In the physical case of A^ — 2, the above first order 
estimates for the critical exponents {v w 1.14, rj ra 0.73) 
compare favorably with the available Monte Carlo results 
{v — 1.00, rj = 0.75) [||, since the actual (inverse) pa- 
rameter of this expansion A^Tr(l(g)l) = 4Af remains fairly 
large even for N — 2. For comparison, the first-order es- 
timate for the anomalous dimension exponent given by 
the e-expansion of Ref. Q is 77 « 4e/7 = 0.57, although 
the agreement can be improved in the higher orders |q]. 

We emphasize that the above discussion pertains to 
the bulk properties of the layered d-wave superconduc- 
tors and is, therefore, unrelated to the previously pro- 
posed scenario of a surface-induced is pairing pj. 

The alternate case of idxy pairing corresponds to a 
different choice of the coupling matrix A = 1 (g) CT3 in 
Eq.(y). This yields a chiral invariant order parameter 

vPal ® 0-3*0 = '4>ia^2'4>ia — i'la'^i^^a which causcs the 
term (pi) to explicitly break the chiral symmetry (m. 



However, contrary to the is order parameter, the idxy 
pairing and the corresponding boson field (j) change their 
signs under parity transformation {x — > —x,y — > y) 
which acts upon the four-component spinors as *„ — > 
CTi (g) CTiVpQ. It is this very property that allows for the 
existence of a zero-field spin or thermal Hall current in 
the idxy-pa,iTed state |^. 

Therefore, by substituting parity for the chiral trans- 
formation (H), one finds yet another Z2 symmetry which 
gets spontaneously broken upon the onset of the idxy 
order described by the same gap equation (||). The re- 
sulting critical behavior then appears to be identical to 
that found in the case of the is pairing. It will be, how- 
ever, different from the previously discussed second order 
d —f d + id' transition induced by magnetic impurities |^ 
as well as the first order one driven by external mag- 
netic field ||l^ (both these scenarios have recently been 
revisited and critically assessed in Ref. [|lj). 

After having fully identified the nature of the sec- 
ond pairing transition, we now proceed with comput- 
ing the experimentally measurable quasiparticle damp- 
ing. Whether attainable by varying a single physical pa- 
rameter such as doping or not, a zero temperature quan- 
tum critical point affects the finite temperature dynamics 
of the system in a whole domain of the g — T plane where 
the correlation length L^ exceeds the thermal one ex \/T 
(the dynamical exponent equals unity due to the Lorentz 
invariance) . 

In this quantum-critical regime, the renormalization 
cut-off is set by the temperature, and, therefore, the small 
energy- and momentum-dependencies of the fermion and 
boson propagators are no longer governed by the anoma- 
lous operator dimensions from Eq.(|8|). Instead, the latter 
appear in the residues oc T^([*^'^l^~of these (now, pole- 
like) propagators. 

After factoring out the fermion wave function renor- 
malization factor, the damping of the nodal quasipar- 
ticles with p = is given by the expression: I](e) = 
Im{Tr[ro(Zp(v[/*))-i]}. 

The universal function E/T — <^{x^y) incorporates, 
alongside with the low- versus high-energy asymptotics 
for x — e/T < 1 (> 1), the crossovers from the quantum- 
critical to the other, renormalized classical and quantum 
disordered, regimes for y = ^{\g — gc\lgcY /T < 1 (> 1), 
respectively. 

The non-linear equation determining the function 
$(2:, y) results from including S(e) into the fermion prop- 
agator in Eq.(6) 
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and, correspondingly, into Eq.(7) for the spatially 
isotropic, albeit manifestly non-Lorentz- invariant, finite 



temperature polarization operator 

n(.,<,)=,"»/|i/ii^{.a.h 
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(e + ^S(£))(£ + UJ + iT,{e + u)) - p(p + q) 
[(e + iE(e))2 - p2][(e + c^ + il](e + tu)y - (p + q)2] 

(10) 

Notably, in the quantum-critical regime (y < 1 ) 
the energy-independent renormalization factors drop out 
from Eqs.(p|,p^), thanks to the underlying Ward identities 
of the HY model g. 

Both analytical and numerical, analyses of the coupled 
equations (§,|l^ show that the fermion damping behaves 
as E ex T and oc e at high and low temperatures, respec- 
tively, in agreement with experiment ||l| and the results 
ofRef. g. 

Below the crossover ( y > 1 ) to the quantum disor- 
dered (weak coupling, low-T) regime, the infrared cut- 
off is provided by the (inverse) correlation length L7^, 
and the damping of the Lorentz-invariant nodal quasi- 
particles appears to retain its energy dependence even for 
e < T, unlike the case of fermions with extended Fermi 
surface. In the energy intervals x > 1, l/j/'* < x < 1, 
and X < 1/y'', the self-consistent solution of Eqs.(9,10) 
behaves as $(x,y) ex x^/y"^, ex x^^^/y"^, and oc 1/y^, 
respectively. For e ~ T it reproduces the perturba- 
tive second order result S ex T'^ (see |^) valid for a 
generic short-ranged coupling, while for e — > we obtain 
S ex T^, one possible experimental implication being a 
significant narrowing of the nodal ARPES linewidth be- 
low the crossover. 

Likewise, in the renormalized classical (strong cou- 
pling, low-T) regime one can expect a suppression of the 
nodal quasiparticle density of states that can manifest it- 
self in tunneling and specific heat data as well as a sharp 
decrease of both, thermal and optical, conductivities ac- 
companied by an increasing thermal Hall angle fl^ . 

Furthermore, according to experiment, the linear tem- 
perature/energy dependence of the quasiparticle damp- 
ing extends well into the normal phase [l|. This might 
be suggestive of a possibility of applying the HY model 
to the pseudogap state, as the above discussion relies on 
the presence of the dj,2_j^2 -symmetrical (pseudo)gap in 
the local quasiparticle spectrum rather than the onset 
of the superconducting coherence across the entire sys- 
tem. To this end, it should be interesting to contrast the 
quantum critical scenario described in this paper against 
the results of the previous approaches to the pseudogap 
problem which have focussed on the couplings between 
the nodal Dirac fermions and magnetic order parameters 
[ fi3[ as well as classical pj] and quantum |l5[ supercon- 
ducting phase fluctuations. 

To summarize, we carried out a non-perturbative anal- 
ysis of the secondary pairing transition in planar d-wave 



superconductors. By identifying the transition in ques- 
tion with the HY model of the 2 -I- 1-dimensional Dirac 
fermions we succeeded in finding all the critical exponents 
and determining the behavior of the quasiparticle prop- 
agators. At e > T the quantum-critical fermion propa- 
gator was found to have a numerically small anomalous 
dimension, while in the opposite limit it exhibits linear 
damping, as observed in experiment. 

Should a secondary pairing occur in the bulk d-wave 
superconductors upon tuning the phenomenological pa- 
rameter (7 in a (yet, unspecified) way, it can manifest itself 
as a marked change in quasiparticle spectral and trans- 
port properties even in the absence of boundaries (cf. [[7|), 
magnetic impurities (cf. ||9|), and/or external magnetic 
field (cf. |1^]). Also, the two-dimensional spatial parity's 
remaining intact or being spontaneously broken may al- 
low for an experimental discrimination between the two 
different kinds {is versus id^y) of incipient ordering with 
broken time-reversal symmetry. 
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